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Abstract 

It has been previously shown that calculation of renormalization group (RG) 
functions of scalar 0^ theory is reduced to thermodynamic properties of 
the Ising model. Using high temperature expansions for the latter, RG 
functions of the four- dimensional theory can be calculated for arbitrary 
coupling constant (?, with an accuracy of 10~^ for the /3-function and with 
an accuracy of 10^^ - 10^^ for anomalous dimensions. The expansions of 
the RG functions up to the 13th order in g'^^"^ have been obtained. 



1. INTRODUCTION 

As was recently shown in [1, 2], the Gell-Mann - Low function [i{g) and anomalous 
dimensions of the 0^ theory can be expressed in terms of the functional integrals, providing 
the representation 

g = F{go,mo,A) , /3{g) = Fi{go,mo, A) , (1) 

where go and mo are the bare charge and mass, respectively; A is the momentum cutoff 
parameter, and g is the renormalized charge. Large g values are reached only near a zero 
of one of the functional integrals, where the right-hand sides of Eqs.l are significantly 
simplified and the parametric representation is resolved in the explicit form. As a result, 
asymptotic expressions for the (3 function and anomalous dimensions are obtained. A 
similar approach can also be implemented in QED [3]. 

Parametric representation (1) has the following general property. If 5^0 is expressed in 
terms of g using the first of Eqs.l and the resulting expression is substituted into the second 
equation, the dependence on mo and A disappears according to the general theorems [4] , so 
that the (3 function depends only on g. However, this property is not automatically satisfied 
in applied calculations. The reason is that the general theorems imply the continual limit 
A — > 00, which physically means the condition 

m < A or ^ > a , (2) 
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where m is the renormahzed mass, ^ is the correlation radius, and a — is the constant 
of a lattice at which the functional integral is defined. Under condition (2) in the region 
of large go values, the functional integrals of the 0^ theory are reduced to Ising sums; as a 
result, Eqs.l have the form 

g = F{k) , /3{g) = F,{k) (3) 

where k has the meaning of inverse temperature in the Ising model, and it is obvious that 
the (3 function depends only on g. Condition (2) formally corresponds to the inequality 
—gQ^ml/A'^ ^ 1, but the reducing to the Ising model is really possible under the weaker 
condition 

^0 > 1 , -go^^VjA'^ > 1 , -gQ^ml/K^ ^k — arbitrary . (4) 

For this reason, parametric representation (3) remains valid in the region of small k, where 
values of g are large and gradient expansions arc applicable. At first glance, the condition 
go ^ 1 corresponds to the strong coupling regime and parametric representation (3) is 
limited by only this condition. However, there is another view on this situation. Let us 
strengthen conditions (4) by passing to the limit 

5^0 — > oo , —gQ^^'^ml/ A'^ ^ oo , — gQ^rn^/ ~ k — const (5) 

In this case, the transition from Eqs.l to Eqs.3 is valid without any approximations and 
conserves strict equivalence with the initial 0"^ theory under a certain choice of its bare 
parameters; the last property ensures the conservation of the form of the Lagrangian under 
renormalizations. The passage to the limit go ^ oo does not mean the same passage for 
the renormahzed charge g\ in fact, according to gradient expansions, g varies from infinity 
to about unity when k varies from zero to about unity. Since parametric representation 
(3) is exact and specifies the (3 function in the interval 1 < (7 < oo, it can be analytically 
continued and treated as a definition of (3{g) at arbitrary g values. However, there is a 
question: Does this definition provide correct results in the weak-coupling region? 

An answer to this question can be obtained using high-temperature scries [5]. Such 
series are traditionally constructed for quantities X21 f^2, ^4 (see Section 2), which com- 
pletely specify the right-hand sides of Eqs.3. High-temperature expansions are formally 
applicable for small k, but their comparatively large length (up to 30 terms in some cases) 
allows a successful analysis of the vicinity of the phase transition point Kc and leads to 
the results consistent with other methods. Consequently, good approximations for the in- 
dicated quantities can be obtained throughout the interval < k < Kc- The substitution 
of such results into the right-hand sides of Eqs.3 makes it possible to determine the renor- 
malization group functions in the interval g* < g < 00, where g* is the fixed point of the 
renormalization group. In the four-dimensional case, g* — and the mentioned procedure 
completely determines the renormalization group functions. In many works [6 — 16], the 
high-temperature series were used to test logarithmic corrections to scaling [17]. Already 
those works provide the positive answer to the above question: parametric representation 
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(3) gives correct results in the weak-coupling region. Therefore, we can concentrate our 
efforts on constructing the rcnormalization group functions of the four-dimensional 0^ the- 
ory for arbitrary g values. This can be done with an accuracy of 10~^ for the f3 function 
and with a slightly lower accuracy for anomalous dimensions. 

The determination of calculated rcnormalization group functions implies the use of a 
lattice regularization different from the usual Pauli-Villars regularization scheme, isotropic 
cutoff in the momentum space, dimensional regularization, etc. However, the /3 function in 
the used scheme is determined in terms of the observed charge and mass [1, 2] and should 
be independent of the cut-off procedure. Such a dependence is possible for anomalous 
dimensions, because they are determined in terms of the unobservable Z factors. In any 
case, the distinction of this way of regularization from the usual procedures is no more than 
difference between the latter procedures 



2. INITIAL RELATIONS 

Let us consider the n component 0^ theory with the action 

[ a=l a=l \a=l J ) 

uo = S'oA^ , e = 4 - d , (6) 

where go and mo are the bare charge and mass, respectively; d is the dimensionality of 
space; and A is the momentum cutoff parameter. The most general functional integral of 
this theory contains M multipliers of the field in the pre-exponential factor, 

^I^.^mI^I'-'-'^m) = j D(p(paAxi)^a2{x2) ■■■(PaM{xM)exp{-S{^}) , (7) 

and will be denoted as KM{Pi} after the transition to the momentum representation and 
the separation of S factors, 

(8) 

where Iai...aM the sum of terms 50,^02^0304 • • • with all possible pairings, and J\f is the 
number of sites of the lattice on which the functional integral is defined. The integrals 
KmIPi} are usually estimated at zero momenta and only one integral i^2{p} is required for 
small p values, 

K2{p)^K2-k2p^ + ... (9) 

Below, the case with d — A and n = 1 is considered, but the general formulas are written 
for arbitrary d and n values. 
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The below consideration concerns the renormahzation group functions P{g), r]{g), and 
7/2(5') entering into the Callan-Symanzik equation [1] 



r(i>^) = o, 



(10) 



for the vertex r^^'^^ with external lines of the field cj) and L external interaction lines. 
The expression of these functions in terms of the functional integrals leads to the parametric 
representation [2] 



^^^^ \kJ Kl \ {InK./K^y 



(11) 
(12) 



vig) = 2 



{\nK2/Koy + (In K^/K^y 

{XnK^/k^y 



(13) 



(14) 



{InK^/K^yilnK./K^y 

where primes stand for derivatives with respect to rn^. Under condition (4), the functional 
integral of the scalar theory can be written in the form 

Z^jxi} = (2k)'^ j dip^ v^xi • • • V^XM exp |-k XI -'x-x'V'xV'x' | H '^(v'x - 1) (15) 

and is transformed to an Ising sum over the values ipy^ = ±1. The quantities studied in 
high-temperature expansions are introduced as 



(16) 



x,y,z 



(where superscript c marks the connected diagrams) and coincides up to factors with the 
ratios K^jK^^ and K^jK^ of the functional integrals introduced above; more pre- 

cisely, 

^ = 2rf— = -/o , 



K2 

— = 2kX2 = Kf2{l^) , 
K/^Kq _ 1 X4 _ r ( \ 



(17) 
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where the introduced functions fi{n) will be used below. It was taken into account that 
there is no zeroth term in the expansion of in (see Eq.20 below), so that all functions 
/o(k), /2(/t), and /4(k) are regular and their expansions begin with the zeroth term. The 
substitution of Eqs. 17 into Eqs. 11-14 gives 



g l-K[ln/o(^)]' 



via) = -2k 

V2{9) = - 



-K[\nfo{K)]'' 



{l-K[\nfo{K)]') {l + K[\nMK)]') 

It is easy to obtain the strong coupling behavior for renormalization group functions 
taking limit k — ?• 

m = dg, Vi9) = 0, V2i9) = -^ (g-^oo). (19) 

For a simple hypercubic lattice with the interaction between the nearest neighbors, the first 
terms of the expansion of functions (16) for d = 4 and n = 1 have the form [12] 

= 1 + 16k + 224^2 + . . . 

/X2 = 16k + 512^2 + 33920/3k^ + . . . (20) 

X4 = -2 - 128k - 4672^2 - . . . 

The substitution into Eqs. 18 makes it possible to obtain the expansion of the renormaliza- 
tion group functions in g~'^^^ and, in particular, a more accurate asymptotic expression for 

vig) 

. ^ 16 1 

9 g 

The universality of this asymptotics has not been tested and, strictly speaking, it refers to 
the indicated model. Below, 14 terms of expansion (20) presented for n = 1 in tables 5, 8, 
and 11 of the paper [18] are used. 

3. VICINITY OF THE PHASE TRANSITION 
3.1. General Strategy 



5 



The foundation of the apphcation of high-temperature expansions for investigating the 
critical behavior is as follows. Let a certain quantity F{n) has a power-law behavior near 
the transition point Kc = '^/Tc 

F (X {T - T^Y^ (X - k)-^ . (22) 

In this case, the convergence radius of the expansion in n is limited by the quantity Kc. 
In actual the nearest singularity to the coordinate origin; this circumstance 

facilitates its analysis. It is easily seen that the nearest singularity for the logarithmic 
derivative 

(lnF)' = :^~-^ (23) 

is a simple pole with a residue —A and can be investigated using the Pade- approximation. 
The Pade-approximant [M/N] is defined as the ratio of the polynomials of the degrees M 
and iV, 

QNi^-) I + qiK + . . . + q^K^ 

whose coefficients are chosen such that the first M + N + 1 coefficients of the expansion of 
(InF)' in k are reproduced. It is known that Pade-approximants successfully predict the 
nearest singularities of the approximated function if these singularities are simple poles [5, 
19]. Diagonal (M = A^) or quasidiagonal (M ^ N) approximants are usually used for which 
convergence to the corresponding function is proved under the most general assumptions. 
The use of this strategy in the four- dimensional case is complicated by the existence of 
logarithmic corrections to scaling [T71 H] : 

X2 ~ r-i|lnr|P, 

where r ~ (kc — k) is the distance to the transition and the exponent p is determined by 
the first terms of the expansion of the renormalization group functions, 

m = f32g' + l^sg' + ..., 

r]{g) = 629' + 6,g^ + . . . , (26) 
V2{9) = Ci9 + (29^ + ■■■ , 

where 

R - Q R - 02 9n+42 _ 2^+2 /- _ q ^+2 .^^n 

P2 — ^4 ; P3 — —^4 Z , O2 — O4 — — , C,i — — — — [Z() 
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and ^4 = I/Stt^. According to (25) we have for functions fi 

/o~r|lnr|-^ f, ^ T-'\\nT\P , /4 ~ r'^l Inr^f-i . (28) 
The behavior of the charge g is given by the expression 

9 = 7^, Co = 2/ (32 (r^O), (29) 
I m r| 

where the coefficient of the logarithmic factor is universaL When Eqs. (28) and (29) 
are vahd, parametric representation (18) automatically ensures the results (i{g) = (329"^, 
ri{g) = ■ g, and ri2{g) = dg, the correct behavior of the renormalization group functions 
at small g values. 

The objective test of Eqs. 25) for lattice models were performed in many works [6-15]. 
In particular, it was convincingly shown in [6, 7] that high-temperature series for the Ising 
model allow reliable prediction of the exponent p. Expression (29) was confirmed with 
a satisfactorily accuracy in [7, 9]. Already these results provide the positive answer to 
the question formulated in the Introduction: parametric representation (18) gives correct 
results for the renormalization group functions in the weak-coupling region. 

3.2. Zeroth Approximation 

The Pade-analysis of Eqs. 28 is performed by the successive approximation method. 
In the zeroth approximation, the logarithmic factors are ignored and the functions fi are 
processed under the assumption of their power-law dependence on k. The results of such 
an analysis presented in Table 1 show a significant difference of the obtained exponents 
from the exact values (see Eqs. 28) and provide a rough estimate of the critical point 

= 0.07476 ^ 0.07490 

A more accurate estimate of Hc can be obtained taking into account that the ratio X4/X2 ~ 
74/2 in the scalar case (when p = 1/3) behaves as and contains no logarithms [6]. As 
is seen in Table 2, the Pade-analysis of this quantity provides the exponent really close to 
the exact value and the corresponding estimate of Kc 

= 0.07481 ^ 0.07487 (30) 

is almost final and will be only slightly refined below. The central value of interval (30) 
almost coincides with the result Kc = 0.074834(15) obtained in [6] with a more sophisticated 
processing. 

3.3. First Approximation 
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In this approximation, the following representation is used: 

/o = /o|lnr|-^ /2 = /2|lnrr, f, = hl^nrl'"-' (31) 

and the Pade-analysis is applied to the functions /j. Since the relation r = A[kc — i^) 
includes the nonuniversal factor A, it can be accepted that 

I In r| = - ln(l - k/kJ , (32) 

Table 1. Position of the pole corresponding to the critical point Kc and residue at 
it (in parentheses) for the Pade approximant [N/N] of functions [ln/j(K)]'.0 





[In/o 




[ln/2(«:)]' 


[ln/4 






2 


0.07519 


(1.130) 


0.07510 (-1.113) 


0.07442 ( 


-1.^ 


B32) 


3 


0.07521 ( 


1.131)* 


0.07543 (-1.085) 


0.07419 ( 


-1.^ 


S14) 


4 


0.07502 


(1.116) 


0.07497 (-1.101) 


0.07476 ( 


-lA 


B79) 


5 


0.07480 


(1.063) 


0.07513 (-1.103) 


0.07477 ( 


-l.l 


B81) 


6 


0.07486 


(1.082) 


0.07490 (-1.088) 


0.07476 ( 


-lA 


B79) 



Table 2. Position of the pole corresponding to the critical point Kc and residue at 
it for the indicated Pade-approximants of function [ln/2/4]'. 





[A^ + 1/N] 


[N/N] 


[N/N + 1] 


2 


0.07418 (-2.871) 


0.07461 (-2.936) 


0.07558 (-2.963) 


3 


0.07488 (-2.993) 


0.07450 (-2.923) 


0.07465 (-2.946) 


4 


0.07486 (-2.988) 


0.07485 (-2.986) 


0.07486 (-2.988) 


5 


0.07487 (-2.989) 


0.07486 (-2.987)* 


0.07491 (-2.998)* 


6 


0.07481 (-2.970) 


0.07484 (-2.983) 


0.07483 (-2.978) 



where the free parameter Aq and trial value for the critical point are used to accurately 
fit the exponent and to obtain a self-consistent result for ■ According to Table 3, such a 

^ The asterisk in Tables 1-4 marks defective approximants. A "defect" in the Pade-analysis is the ap- 
pearance of a pair of a pole and a root close to each other; as a result, the corresponding Pade-approximant 
is reduced to a lower order approximant. The defectiveness of the approximant can lead to loose of the 
accuracy and is a reason for its discrimination. 
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fit is easy and good results for the exponent are obtained in a wide range of the Aq values. 
The optimal Aq lie in the interval 0.13 -j- 0.63 and a new estimate of the critical point 

= 0.07483 ^ 0.07489 

Table 3. Pade analysis of the functions /j, introduced according to Eqs.31. 















Ao 


[ln/o(«:)]' 


[ln/2( 


k)]' 


[ln/4( 






[6/6], = 0.074842 


[6/6], = 


0.074834 


[6/5], = 


0.074890 


2.0 


0.07491 (1.037)* 


0.07486 ( 


-1.023) 


0.07488 (- 


-1.968)* 


1.0 


0.07487 (1.018) 


0.074844 ( 


-1.007) 


0.07493 (- 


-1.989)* 


0.625 




0.074834 (- 


-1.00005) 






0.5 


0.074855 (1.0085) 


0.07482 ( 


-0.996) 


0.07477 ( 


-1.960) 


0.25 


0.074846 (1.0029) 


0.07475 ( 


-0.973) 


0.07488 ( 


-1.996) 


0.2 








0.074890 ( 


-1.9994) 


0.13 


0.0748420 (1.00005) 










0.1 


0.074840 (0.9993) 


0.07483 (- 


-0.990)* 


0.07490 (- 


-2.0044) 


0.06 


0.07487 (1.0033) 


0.07482 (- 


-0.988)* 


0.07491 (- 


-2.0063) 



is only slightly shifted as compared to Eq.30. The results for the constant Cq in Eq.29 are 
shown in Fig. la; they are close to the theoretical value Cq''* = 35.09, but are systematically 
above it. Similar inaccuracies in the determination of Cq were observed in other works. The 
use of constants A, B, and D for a simple hypercubic lattice from Table 5 in [7] gives the 
estimate cq = B/A'^D'^ = 142.8 instead of a theoretical result of 105.2 referring to the used 
normalization. A worse estimate was obtained in [10]; very bad results (discrepancies of 9 
and 18 times) were obtained for other lattices [7]. A satisfactory test of Eq.29 was declared 
in [9], where the tested relation was not Eq.29, but its consequence dg~^/d\n.T — I/cq; 
in this case, the central value cq approximately corresponds to Fig. la and the agreement 
with the theory was achieved at the expense of an increase in the uncertainty of the results 
because of differentiation. 

3.4. Second Approximation 

Expressions (25) and (28) are obtained in the leading logarithmic approximation. In 
the next-to-leadinglogarithmic approximation (see Appendix A), they have the form 

fo = hoT {fsing)~^ , /2 = ^2 {fsingY hging , A = ^4 {fsingY^~^ , (33) 

Here, the functions hi{K) are regular at k — > and singular functions are chosen in the 
form 

fsing{n) = 1 - ^ In T In (1 - ^ In r) , (34) 
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Figure 1: Constant Cq in Eq.29 versus the parameter Aq in the leading logarithmic approx- 
imation (a), and versus g in the next-to-leading logarithmic approximation (b). 
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hsing{l^) = 1 + 7-T , T ^ 1 - k/ Kc (35) 

where 

_ 2/33 Ci _ n2-8n-68 _ 2^2 _ n+2 

{n+sy ' ^"^"(^rPsF- ^^^^ 

The main distinction from Eqs.28 is reduced to the replacement of | In t| by | In r| +s In | In r | 
with a known parameter s; in view of the ambiguity of the normahzation of r, it is necessary 
to consider the combinations A + | lnr| + s\n{B + | lnr|). where the constants A and B 
are different for different functions. Formally, these constants do not affect the character 
of a singularity, but their unsuccessful choice can strongly distort the results. To avoid 
a large number of fitting parameters, fsing{i^) was taken in the functional form following 
from perturbation theory. A reason for such a choice is as follows. The parameter g has 
the sense of the Ginzburg number and determines the size of the critical region, where 
logarithmic corrections are significant. It is of interest to estimate this parameter, because 
the Ginzburg number is often small even in the absence of theoretical reasons for this. The 
function fsing{i^) at small values of g is close to unity almost everywhere, but increases 
sharply near Kc- If the singularity is separated inappropriately, regular functions hi{K) 
in Eqs.33 are rapidly varying near and are poorly reproduced by Pade-approximants. 
However, for small values of g the form of Eq.34 is practically exact, so that functions hi{K) 
are almost constant. For ^ ^ 1, the form of Eq.34 is not exactly correct, but inaccuracy in 
the separation of singularities in this case is not so critical, because the function fsing{i^) is 
a rather slowly varying. 

The universal choice fsingi^^) foi' functions is possible if the 0{g) contributions are 
negligible as compared to unity (see Appendix A), so that the inclusion of factors of the 
hsing{i^) type is strictly speaking beyond of accuracy. However, such factors are sometimes 
of qualitative importance. In Eqs.33 they are taken into account in the minimal manner: 
the product /0/2 in this form has the correct singularity and ensures the correct behavior 
of r]{g) at small g\ similarly, the product 74/2 is incompletely free of logarithms and this 
property makes it possible to slightly correct deviations observed in Table 2. 

Table 4 presents the Pade-analysis of the functions fi introduced by the relations 

/o = /o if sing) ^ , /2 = /2 {fsingY hging , f^ — fi {fsingY^ ^ , (37) 

rather than by Eqs.31; the estimate of the parameter cq in Eq.29 is illustrated in Fig. lb. 
It is easily seen that the actual interval of g vahies is much narrow than that in the leading 
logarithmic approximation (wherethe parameter 1/Aq is similar to g). The optimum values 
for various functions cover the range of 0.85 -j- 1.06, which provides the estimate 

Co = 36.3 ±1.8 (38) 

in good agreement with a theoretical value of 35.09. The exact cq value is realized at 
g 1.02 (see Fig. lb). Finally, Table 4 presents the maximally accurate estimate of the 
critical point 

Kc = 0.074840 0.074867 , (39) 
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which is available with the existing information. The values accepted below are Kc — 
0.074850 from the middle of interval (39) and g — 1.020385, which ensures the exact Cq 
value for the [3/3] approximant. 



Table 4. Pade-analysis of the functions /j introduced according to Eqs.37. 
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[ln/o(«)]' 


[ln/2(«)]' 


[ln/4(«)]' 




[6/6], = 0.074843 


[6/6], = 0.074840 


[6/5], = 0.074867 


0.5 


0.07492 (1.036)* 
0.07488 (1.019)* 


0.07488 (-1.024) 


0.07487 (-1.968)* 
0.07491 (-1.988)* 


0.7 


0.07485 (-1.0096) 


0.85 




0.074840 (-1.0008) 




0.9 


0.07485 (1.0052) 


0.074836 (-0.998) 


0.074877 (-1.994) 


0.99 


0.074843 (1.00005) 






1.0 


0.074842 (0.9995) 


0.07483 (-0.994) 


0.074865 (-1.997) 


1.06 






0.074867 (-2.0001) 


1.2 


0.07482 (0.988) 


0.07476 (-0.976)* 


0.07488 (-2.010) 



4. RESULTS FOR RENORMALIZATION GROUP FUNCTIONS 

The derivatives of singular functions can be written in the form 



] ) [1^ fsinq\ / \ 9 ) [1^ hsing\ ) [ln ^sing] / \ 



where 



ui{t) = - — 1 + 



fsing if sing + QQ) \ 1 " 5' In T / ' 

V2{r) = 77 ■ =r < 1 + 



fsingifsing + QQ) [ 1 - gluT {1 - gin t)^ 

/ 9 ^ 9 _ ](, , s-g 



Jsing fsing + q9 J \ 1 " 5' In 
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(40) 



fsing V 1 - 5' In T y ' 

Mr)=^ |l+ ^^- (^ '-^1 ^2 |-^l(^)^ (41) 
fsing [ l-glnr {1 - glnry } 

( N <lf (. . -^9 \ 



Taking into account (40), the substitution of Eqs.33 into Eqs.18 provides the parametric 
representation for the renormahzation group functions in the form 



H(k) 



H(k) 



1 



k2 f . ' 

^ J smg 

P{g) 2Kcr {2- K [In h^hl]') + 2kui 



vig) 



KcT — K [Inho]') + K{l+pUi) 

—2kcTK [In /?.o^2]' — 2kvi 
KcT — K [In ho]') + k{1+pui) 



(42) 



[KcTf {l-K^ [In h2]") + {KcT{1 + K [ln/l2]0 + K{l+pUi+Vi)f - k'^{1+PU2 + V2) 

[kcT (1-k [In/io]') + k(1+pmi)} [k^t (1 + k [ln/i2]') + K{l+pui+Vi)} 



Asymptotic expressions (19) are obtained at k — irrespective of the form of regular 
functions, whereas at r — )■ we have the results 

2^ m 2g 2{s-p)f 2qf 2pg 

3=ir-r — > — + — f2 ' v{g)=-P2—^ V2{g)=—. — , 

P2 J sing y J sing J sing J sing J sing 

(43) 

which reproduce the first two terms of the expansion for P{g) and the first terms of the ex- 
pansions for ri(g) and 'r]2{g) in Eqs.26.[l When the terms with r are neglected, Eqs.42 provide 
the regular expansions of the renormahzation group functions in g (certainly without the 
reproduction of correct coefficients), whereas the terms with r provide the exp{— const /g) 
singularity, which should exist owing to the factorial divergence of the perturbation series 
[20, 21]. Thus, the parametric representation is rather "intelligent" and ensures the correct 
analytical properties at g ^ 0. 

The accuracy of the entire construction is determined by the accuracy of the determi- 
nation of the regular functions hi{K,). The expansions of these functions in k are obtained 
from Eqs.33 and are used to construct the Pade-approximants, which are regular in the 
interval (0, Kc), because all singularities have been separated. The obtained regular func- 
tions are shown in Fig. 2. For the functions H{k,) and [In /io(/t)]', all approximants provide 
almost coinciding results; small distinctions are visible for the function [ln/i4(K)]' near Kc 
(see Fig. 2). The situation is less satisfactory with the function [ln/i2('^)]' for which an 
increase in the order of the Pade-approximation leads to an increase in the deviations from 
the regular behavior predicted by lower approximants. It is unclear whether the sequence 
of approximants converged sufficiently or such deviations will further increase. Moreover, 

^ Note that the coefBcients /32, /Ss, S2, Ci exhaust invariant (scheme- independent) information on the 
renormahzation group functions and a further refinement of the procedure (the construction of the next- 
to-next-to-leading logarithmic approximation, etc.) requires the calculation of the subsequent coefficients 
for the corresponding lattice regularization. 
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Figure 2: Regular functions H{k) and [ln/ij(K)]' obtained in the Fade- approximation. 

these deviations can be artifact due to an incompletely consistent separation of singulari- 
ties leading to a residual singularity in the function [ln/t2(«^)]' (in the used approximation), 
which affects higher approximants. In the latter case, the behavior predicted by the [3/3], 
[2/3], and [3/2] approximantscan be more authentic. Fortunately, this dilemma can be 
resolved using the strong-coupling expansions (see Section 5), which certainly indicate that 
the use of higher Fade-approximants is correct and the results obtained in this case are 
satisfactory. Appendix B presents the parameters of the approximants used for H{k,) and 
[ln/ij(K)]', which allow the application of parametric representation (42). 

To represent the results, it is convenient to use the so called "natural normalization" of 
the charge, which is obtained by the change g — > (167r^/3)(7 and corresponds to the repre- 
sentation of the interaction termEl in the form (167r^/4!)(7o0^; in this case, the parameter a 
in the Lipatov asymptotic form ca^r{N + b) [20, 21] is unity and the nearest singularity 
in the Borel plane lies at the unit distance from the coordinate origin [21]; this property 
defines functions varying at an approximately unit scale. The solid lines in Fig. 3 are the 
resulting renormalization group functions, whereas the dashed lines are the strong- and 
weak-coupling asymptotic behaviors. The approach to the strong-coupling asymptotics is 

^ The traditional representation g^cfA /9> in the n-component case is motivated by the fact that the vertex 
^a^^7(S ~ gla/B-yS in thc lowcst Order is golapjs^ which ensures the relation g — go in the limit go — > 0. In 
the scalar case, the tensor lap-yS is reduced to three and the interaction is represented as 5o0^/4!- This 
motivation logical at first glance is in fact illusory, because the bare charge has no physical sense. 
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strongly prolonged in agreement with the results reported in [22]. However, the prolonga- 
tion of the one-loop behavior of the /3 function pointed out in that work is not confirmed: it 
appears to be an artifact, conditioned by essential exceeding of the limiting value of f3{g)/g 
obtained in [22] in comparison with Fig. 3 [1] . 

To illustrate the accuracy of the construction, the dotted lines show the results obtained 
if the functions are changed to constants; in this case, the results contain no infor- 

mation on these functions, because [ln/ij(K)]' = and a constant value of H{k,) is fixed by 
Eq.29. It is easy to see that an accuracy of about 1% for (3{q)/g and ri2{g) is reached even 
in the complete absence of information on regular functions. □ 

The real uncertainty of the construction is about two orders of magnitude smaller than 
the difference between the solid and dotted lines, because the regular functions (see Fig. 2) 
are specified better than 1% except for the region k > O.Skc, where the error for the function 
[In /;,2(/t)]' = can reach 10%. However, this region corresponds to g < 0.5 (see Fig. 4), 
where the effect of regular functions is insignificant. 

Table 5. Coefficients of the expansions in g~'^^'^ for the functions I3{g)/g, r](g) and 

V2{g)- 



N 




vig) 


V2{g) 





4.0000000000000 


0.0000000000000 


-4.0000000000000 


1 


-26.127890589687 


0.0000000000000 


26.127890589687 


2 


106.66666666666 


1.7777777777777 


-60.444444444444 


3 


-557.39499924665 


-11.612395817638 


81.286770723472 


4 


3214.2222222221 


29.708641975308 


-44.879012345695 


5 


-16396.702894504 


22.708685154477 


-1208.7213779957 


6 


67356.444444432 


-961.13125612398 


9071.1992161454 


7 


-139720.34647768 


7188.4949076856 


-49662.878604241 


8 


-717634.37037244 


-27680.892323840 


197619.39191503 


9 


9878174.8209247 


-7609.7703277375 


-226822.08364126 


10 


-59767955.489704 


938372.27840847 


-3873286.8465521 


11 


186179701.36334 


-7226487.6363735 


41826925.334797 


12 


355069103.58896 


27981910.625966 


-249549251.38460 


13 


-8851453360.7421 


7407298.5714308 


794136522.54618 



The reason is that the terms [In/i^]' in Eqs.42 has the factor ) , which is smah both for 

K — and for k — > k^] this factor in the middle of the interval k = Kc/2 is equal to k^/^i whereas the other 
terms are on the order of Kc- In view of Kc ~ 1/15, the effect of regular functions on l3{g)/g and 772(5) is 
about 1%. The situation for 77(17) is different in view of the absentof the KcT term in the numerator. 
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Figure 3: Solid lines are the renormalization group functions. The dashed lines arc the 
strong- and weak-coupling asymptotic behaviors. The dotted lines are the results obtained 
under the assumption of the constancy of regular functions hi{K,) under which Eqs.42 con- 
tain no information on them. 
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5. STRONG-COUPLING EXPANSIONS 



Expanding the right-hand sides of Eqs.18 in k and expressing k in terms of g, it is easy 
to verify that the functions l3{g)/g, r]{g), r]2{g) are expanded in g~'^^'^ as 

M=f (44, 

3 N=0 

The expansion coefficients up to iV = 13 recalculated from high-temperature series are 
given in Table 5. El 

It is easy to verify that the ratios B^j^i/Bn are the same order of magnitude for all N, 
indicating the finite convergence radius. The Pade-analysis of series (44) reveals poles in 
the region ~ 0.1; these poles for most approximants do not lie on positive semiaxis 

in agreement with regularity of the renormalization group functions. To obtain the correct 
power-law behavior in the limit (7 — )■ 0, it is necessary to use the [N/N + 2] approximants 
for f3{g)/g and 772(5') and the [N/N + 4] approximants for ri{g). Such a procedure predicts 
62 with an accuracy of about 20%, whereas (32 and Ci are estimated only by the order of 
magnitude. For this reason, the summation of series (44) in the region of small g gives less 
accurate results than the procedure described above. 

All approximants provide almost coinciding results in the region of large g ; this coin- 
cidence holds to g = 0.5 with anaccuracy of about 1%. Such estimates for the functions 
P{g)/g and 772(5') are in agreement with the more accurate results obtained above. The es- 
timates for the function T]{g) certainly indicate that the highest order approximants should 
be used for [ln/i2]' and the results are confirmed at a level of about 1%. Series (44) can 
apparently be used more efficiently, but analysis of this possibility is beyond the scope of 
this work. 



6. DISCUSSION OF THE RESULTS 



The resulting /3 function is non-alternating and has the asymptotic behavior P{g) = Ag 
in the limit 5 — )■ 00. According to the classification proposed by Bogoliubov and Shirkov 
[23] (see discussion in [1]), this means the possibility of the construction of a continual 
theory with a finite interaction at large distances. The last conclusion contradicts the 
widespread opinion that the 0^ theory is "trivial" [24-28]. As was discussed in [1, 30], two 
definitions — Wilson triviality [24] and mathematical triviality [25, 26] — were confused in 
the literature . The first triviality is firmly established (it corresponds to positivity of the 
P function), whereas pieces of evidence in favor of the second triviality are scarce [27] and 
allow another interpretation [1-30] . According to above analysis, we have no contradictions 
in the properties of the lattice 0^ theory with the works cited in [1, 30]. However, there is 

^ Fourteen digits output by a computer are formally presented. The accuracy decreases beginning with 
N — 3 and the last four digits are unreliable at iV = 13. 
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a conceptual contradiction which we want to stress: it concerns the role and significance of 
the lattice theory. 

The usual point of view implies that the lattice 04 theory provides a reasonable ap- 
proximation for the actual field theory. This interpretation provides the natural condition 
^ ^ a, according to which many sites of the lattice should be at the characteristic variation 
scale of the field (f){x). This condition can be hberalized to a or strengthen to ^/a — > oo. 
In the former case the resriction g ^1 for renormalized charge is obtained (for the natural 
normalization) [28], while g = in the latter case (corresponding to the phase transition 
point). Thus, the usual statements are obtained: the theory is trivial in the continual limit 
(A/m — )■ oo), whereas in the presence of a cutoff the interaction is limited from above and 
cannot be strong. The latter circumstance is used to obtain an upper bound for the mass 
of the Higgs boson [28, 29]. 

Our position is that the lattice theory should not be considered as any approximation 
to the actual field theory (although this is possible at go <^ 1). The continual theory 
fundamentally involves no lattice; a lattice appears only in the bare theory, which is an 
auxiliary construction and is completely eliminated later. The bare theory has no physical 
sense and should not satisfy any physical requirements. Without restriction ^ ^ a, the 
renormalized charge can have any value (see Fig. 4). The proposed concept is completely 
consistent with the "rules" accepted in mathematical works [25, 26] according to which the 
continual limit a — )■ is taken at arbitrarily chosen dependences go{(i) and mo (a); in this 
paper, they are taken under conditions (5). 

The only alternative for the perturbative approach is that all quantities referring to the 
continual theory are expressed in terms of functional integrals. These integrals depend on 
go, mo, and A and, with dimensionality taken into account, we have for the charge, mass, 
and other physical quantities Ai (observables, renormalization group functions, etc.) 

g^Fg {go, mo/A) , m ^ KFg {go, mo/A) , = A^F^ (^o, mo/ A) , (45) 

where di is the physical dimension of the quantity A^. According to Eqs.45, the real 
designation of the bare theory is to ensure the representation of the physical quantities 
in a parametric form. The relations between g. m. and Ai are of physical interest; the 
parametric representation is of no deep sense in view of its ambiguity: it can be written in 
various forms by changing go and mo/ A to any other pair of variables. For this reason, 
an attempt to give the physical sense to the bare theory faces the question: Why one of 
numerous parametrizations is of particular significance? 

Excluding go and mo/A in favor of g and m/A, it is possible to arrive at the relation 

Ai = m''^Fi{g,m/A) . (46) 

In the general case, the exclusion of the dependence on A requires the passage to the limit 
m/A — )■ 0, which corresponds to the critical point and returns us to the "zero charge" 
situation. However, the central point is that the general-position situation does not occur 
in Eq.46: after the transformation to the Ising model (vahd under conditions (5)), all 
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Figure 4: Renormalized charge g (a) and mass m (b) versus k/kc- 

functions in Eqs.45 depend on the single parameter k; as a result, the dependence on m/A 
is completely absent[§ in Eq. (46) 

A, = m'^^F, (g) . (47) 

The renormalization program is thereby completed and no additional passages to limits 
are required. This means that (a) the lattice can be retained in the bare theory (as a 
convenient technical tool for the representation of functional integrals) and (b) the relation 
between m and A can be assumed to be arbitrary, which ensures the attainability of any 
value of g (see Fig. 4). 

We consider the above procedure as a real scheme for constructing the continual 0^ 
theory with a finite interaction. In fact, dependence of g and m on bare parameters (Fig. 4), 
as well as the results for the renormalization group functions (Fig. 3), have been obtained 
in the present paper. 

APPENDIX A. Next-to- Leading Logarithmic Approximation 

^ This is not surprising, because the passage to the continual hmit was performed in the process of 
the transformation to the Ising model [2], which was required by the needs of renormalized, but not bare 
theory. 
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The basic formulas referring to the next-to-leading logarithmic approximation under- 
lying representation (33) will be given below. The starting point is the Callan-Symanzik 
equation in the cutoff schemed 



+ p(^o) ^ 7(^0) 



9 In A dgQ 
where the function F satisfies the logarithmic expansion 



F{go,A/m) = 0, (Al) 



/ A \ 

F{go,A/m)=Y.9oT.A§lln-) . (A2) 

^r— n i^— n ^ Til/ 



N=0 K=0 

The substitution of (A. 2) to (A.l) taking into account the expansions 



M=2 M=l 



yields the system of recurrence relations for the coefficients A^: 

N-K+l 

-KA^= E [(3m+i{N-M)-^m]A^zIi, K=1,2,...,N (A3) 

M=l 

In particular, for K close to N 

-NA^ = [P2iN~l)~^,]A^Zl (A4) 

-(AT - 1)A^-' = [(3,iN - 1) - 71] A^zl + MN - 2) - 72] A^zl , 

-{N - 2)A%~^ = [(32iN - 1) - 71] A^zl + mN - 2) - 72] A^zl + - 3) - 73] A^zl , 

etc. The first equation in (A. 4) is solved immediately; after that, the next equations can 
be solved one-by-one using the method of variation of constants. 

Vertex T^^''^\ For this vertex, 'y{go) = i]2{go), all coefficients are nonzero, and Aq = 1; 
the first two equations in Eqs. (A. 4) give 

" ^ TipnN + 1) ' ^ - /32 " /32 ^ ' 



T{l+p)T{N) [ P2^in + p 



(1) 



The substitution of (A. 5) into (A. 2) and the summation of the corresponding series using 
the formulas 



^Its difference from Eq.lO is of no significance at present context, because the first coefficients /32, /Sa, 
§2, Ci ^-I'G independent of the renormahzation scheme. 
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yield 

r(i'2) = |l + 0{go) + ^290 In - + In U + ^29^ In -) 1 (^.7) 
[ m P2 \ mj J 

The 0{9q) terms will be omitted below. 

The renormalized charge g satisfies Eq.(A.l) with 7(5^0) = 0, whereas all coefficients 
in expansion (A. 2) are zero and ^4^ = 1. Similar to Eqs. (A. 5) and(A.7) we have a result 



A^' = {-f32f-' , A^-^ = {-(32f-' (AT _ 1) II i + 0(1) ~ 



(A.8) 



and 



g = goll + P290 In - + In f 1 + P290 In -) 1 (^.9) 
[ m P2 \ mJ J 

which can also be obtained directly from the GellMann - Low equation. 
Renormalized mass.. Neglecting the Z factor, the Ward identity 

p(l,2) ^ p(0,2) ^ _J_ ^ /^ -^QN 

dml drriQ Z 

can be written in the form dml/dvn?' = l/T^^'"^^; the integration with respect to within 
the necessary accuracy is reduced to the multiplication by m^, 

= {ml - m^) ( 1 + ^290 \n- + go^\n(l + P290 In -) 1 , {A.ll) 
[ m P2 \ rnj J 

where ml is the value of m^ corresponding to the transition point. The introduction of the 
dimensionless distance to the transition r oc {m^ — m^) and iterative exclusion of m from 
the righ-thand side give 

= T [1 + ^In l/r + s^ln (1 + ^In 1/r)]-^ , g = ^29ol'^ (A12) 

where s is given in Eq.36. Similarly, (A. 9) reduces to the form 

^ ~ ^l + ^lnl/T + S5ln(l + ^lnl/r) ' ^^'^^^ 

The Z /actor satisfies Eq. (A.l) with 7(5'o) = —v{9o)j while Aq = 1, = A\ = and 
A]^ = for > 2 in expansion (A.2). Similar to (A.8), we have for N >2 

A^' = Al i-P2f-' , A^' = Al i-f32f-' (A^ - 1) I E - + 0(1)1 , (A14) 

L P2 n=2 ^ J 
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and after summation 



4^ (l + ^.,oln A + 90^ In (l + P^g^ln A) ' 

P2 m P2 \ rnj 



(A15) 



Taking into account the relation A\ = —62, expressing m in terms of r and omitting an 
insignificant constant factor, one obtains with the necessary accuracy 



1 + 



2S2 



Pi l + (?lnl/r + s(?ln(l + ^lnl/T) 
The substitution of (A.12), (A. 13), (A. 16) into the relations 



K2 
K2 



m 



K2 



K4K0 



_9_ 



(A16) 



(A17) 



yields Eqs.33 for /^(k). The difference of the Z factor from unity corresponds to the 
corrections of the order go/lnr, which were neglected above, and strictly speaking is beyond 
the accuracy. However, without the inclusion of the Z factor, the product /0/2 would be 
a regular function and, correspondingly, the behavior of ^^(^f) at small g values would be 
incorrect. For this reason, the function hging corresponding to the Z factor isintroduced in 
Eqs.33 by the minimal manner to ensure the correct singularity in /o/2- 

Table 6. Parameters of Pade-approximation (24) of regular functions 

H(k) [\nho{K)]' 



n 


Pn 


qn 


Pn 


qn 





0.166666 


1.000000 


-2.389114 


1.000000 


1 


2.173343 


12.28756 


39.93594 


1.218909 


2 


-8.874246 


-6.056224 


134.2565 


-14.76806 


3 


103.5876 


-124.8396 


-1759.943 


498.1762 


4 








14434.97 


-2468.179 



[\nh2{K)]' 



[In hi{K)]' 



n 


Pn 


(In 


Pn 







2.416517 


1.000000 


5.530725 


1.000000 


1 


-50.63241 


-3.794992 


13.37787 


21.09480 


2 


-345.9676 


-201.7335 


630.6971 


57.28333 


3 


9156.772 


738.3887 


3430.220 


252.1934 


4 


-1285.833 


4787.275 





10511.06 


5 


-267488.9 


-26827.13 








6 


109199.7 


363530.4 
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APPENDIX B. Pade- Approximation of Regular Functions 

Table 6 shows the coefficients p„ and g„ in Eq.24 for the Pade-approximantion of the 
regular functions H{k) and [ln/ij(K)]'; the lowest order approximants having the complete 
accuracy are presented. The singularities were separated with the values Hc = 0.074850 
and g = 1.020385. 
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